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Представлены результаты математического моделирования конвекции вязкой несжимаемой жидкости в прямоугольной поло- 
сти с источниками ввода и отвода массы в сопряженной постановке в рамках приближения Буссинеска. Исследованы режимы 
смешанной конвекции в полости, имеющей две вертикальные и одну горизонтальную стенки конечной толщины, с двумя участ- 
ками ввода и вывода жидкости, и одной свободной поверхностью. Рассмотрена плоская нестационарная задача в рамках моде- 
ли Навье-Стокса для жидкой фазы и теплопроводности для твердой фазы. Получены распределения гидродинамических пара- 
метров и температур при использовании различных граничных условий на внешнем контуре рассматриваемой полости. Устано- 
влено, что даже при умеренном теплоотводе на внешних границах области решения приводит к масштабным изменениям в тем- 
пературных полях жидкости. 


1. Введение 

При решении многих задач теплоэнергетики, 
строительства, химических технологий возникает 
необходимость анализа тепловых режимов объек- 
тов, представляющих собой полость, заполненную 
несжимаемой жидкостью при наличии источников 
ввода и отвода массы, значимых градиентов темпе- 
ратур по всем координатным направлениям и те- 
плообмена по внешнему контуру полости [1, 2]. В 
таких условиях может реализовываться режим сме- 
шанной конвекции, осложненный теплоотводом с 
внешних границ области решения. Во многих 
практически значимых случаях этот теплоотвод 
играет важную роль в формировании теплового ре- 
жима объекта. До настоящего времени моделиро- 
вание таких течений в сопряженной постановке 
[3-6], учитывающей влияние внешних условий на 
характер течения и температурное поле среды, не 
проводилось. Целью данной работы является чи- 
сленное моделирование смешанной конвекции 
вязкой несжимаемой жидкости в прямоугольной 
полости с источниками ввода и вывода массы в 
условиях теплообмена с внешней средой в сопря- 
женной постановке. 

2. Физическая модель 

Рассматривается течение несжимаемой вязкой 
жидкости и теплообмен в полости, имеющей две 
вертикальные, одну горизонтальную стенки конеч- 
ной толщины и одну свободную поверхность, с 
участком ввода и вывода жидкости (рис. 1). 

При постановке задачи принимается, что тем- 
пература вводимой жидкости существенно превы- 
шает начальную температуру среды в полости и 
считается известной. Массовый приход также яв- 
ляется известной величиной. На внешних границах 
области заданы неоднородные граничные условия. 
Неоднородность обусловлена различными тепло- 
физическими характеристиками внешних сред и 
условиями теплообмена на границах области реше- 
ния. Жидкость считается вязкой, теплопроводной, 
режим течения ламинарным. Отток массы с верх- 


ней границы за счет испарения не учитывается. Все 
границы, кроме участков ввода и отвода массы, 
считаются непроницаемыми для жидкости. Тепло- 
физические свойства жидкости и твердой фазы не 
зависят от температуры. 
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Рис. 1 . Область решения задачи: 1) твердая фаза; 2, 3) область 
ввода и вывода жидкости в полость; 4) полость с жид- 
костью; 5) свободная поверхность жидкости 


Задача рассматривается в нестационарной по- 
становке в связи с необходимостью учета способно- 
сти материала твердых стенок аккумулировать тепло 
и значительными перепадами температур между 
вводимой в полость жидкостью и начальной темпе- 
ратурой твердых стенок. В общем случае возможны 
различные варианты расположения отверстий для 
ввода и вывода жидкости. Рассматриваются наибо- 
лее типичные варианты (рис. 1), реализация кото- 
рых приводит к формированию наиболее интерес- 
ных для анализа структур течения жидкости. При 
постановке задачи также предполагается, что кон- 
такт на границах раздела «жидкость - твердая стен- 
ка» идеален, а условия теплообмена на внешних гра- 
ницах области решения не зависят от времени. 

Особое значение в задаче имеют граничные 
условия на внешних границах области решения 
(рис. 1). Эти условия определяют интенсивность 
теплоотвода с внешних границ и могут оказывать 
существенное влияние на режим течения жидко- 
сти. Соответственно, возможно изменение интен- 
сивности теплопереноса в жидкости. Исходя из 
анализа различных вариантов реализации рассма- 
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триваемой схемы на практике можно на правой, 
левой и нижней границах области решения выста- 
влять граничные условия 1-3 или 4-го рода. Но для 
корректного использования условий четвертого 
рода необходимо расширить размеры области ре- 
шения путем ввода дополнительных подобластей 
до значений хи у, при которых фронт прогрева в 
любых режимах не доходит до этих внешних гра- 
ниц. Такой подход существенно усложняет алго- 
ритм решение задачи. Использование условий те- 
плоизоляции на внешних границах области реше- 
ния (рис. 1) не является правомерным, т. к. в реаль- 
ных условиях при различных градиентах темпера- 
тур всегда осуществляется сток тепла во внешнюю 
среду с этих границ. По этим причинам для анали- 
за были выбраны граничные условия первого и 
второго рода на границах внешнего контура. 

Ещё одна особенность рассматриваемой задачи 
заключается в граничных условиях на свободной 
поверхности, как для уравнения энергии, так и для 
уравнений движения. На практике эти условия не 
могут быть определены точно. При постановке за- 
дач математического моделирования решение со- 
пряженной задачи и на этой границе создает прак- 
тически непреодолимые препятствия в связи с 
трудностями описания свободноконвективного те- 
чения газовой (как правило, воздушной) среды 
вблизи свободной поверхности движущейся жид- 
кости. Примеров реализации таких задач пока нет. 
По этим причинам область решения в рассматрива- 
емом случае ограничивалась сверху свободной по- 
верхностью. Предполагалось, что на этой границе 
реализуются условия выполнения граничных усло- 
вий третьего рода для уравнения энергии, а также 
равенство сдвиговых напряжений в газе и жидкости 
на свободной поверхности. В такой интерпретации 
постановка задачи предусматривала возможность 
образования на этой поверхности элементов цирку- 
ляционных течений, обусловленных как естествен- 
ной, так и вынужденной конвекцией. 

При решении свободно-конвективных задач те- 
плообмена чаще всего на внешних границах ис- 
пользуются условия первого рода для уравнения 
энергии [7—10]. Такой подход удобен при числен- 
ной реализации и в целом не выглядит противоре- 
чивым. Но анализ типовых примеров подобных за- 
дач показывает, что условия первого рода в целом 
почти адекватны условиям теплоизоляции. По 
этой причине они не соответствуют реальным ре- 
жимам теплообмена на практике и могут быть ис- 
пользованы только в качестве условий для модель- 
ных задач. Но в этом качестве они являются эл- 
ементами эталонных задач, решения которых игра- 
ют роль своеобразных «измерительных приборов» 
при оценке тех или иных значимых физических 
эффектов на внешних границах областей решения. 
В данной работе оценка масштабности влияния те- 
плоотвода с внешнего контура также проводилась 
путем сравнения с результатами решения такой 
«эталонной» задачи с граничными условиями пер- 


вого рода на вертикальных (за исключением участ- 
ков ввода и вывода жидкости) границах и на ни- 
жней («несвободной») границе. 

Следует также отметить, что моделирование те- 
плоотвода с трех внешних («твердых») границ 
области решения проводилось путем использова- 
ния граничных условий второго рода. При этом ве- 
личина отводимого теплового потока у являлась 
параметром задачи и выбиралась исходя из реаль- 
ного диапазона его возможного изменения. Такой 
прием позволяет не расширять область решения до 
размеров, при которых с одной стороны было бы 
правомерным применение условий четвертого рода 
на границах сопряжения разнородных сред. Но с 
другой стороны такое расширение области анализа 
делало практически невозможным сам численный 
анализ рассматриваемой нелинейной задачи. Вы- 
бор же значений всегда может быть проведен с 
требуемой для практики или для обоснования до- 
стоверности результатов теоретического анализа 
точностью. Так, например, такой выбор может 
быть реализован по итерационной схеме путем 
расчета на каждой последующей итерации величи- 
ны исходя из полученного на предыдущей итера- 
ции градиента температур вблизи границы. Кроме 
этого достаточно точно могут быть определены ве- 
личины исходя из исходных данных каждой за- 
дачи. Последнее создает объективные условия для 
получения достоверных результатов. 


3. Математическая модель и метод решения 


Постановка задачи по пространственным коор- 
динатам является двухмерной и описывается си- 
стемой нестационарных уравнений Навье-Стокса 
для жидкой фазы и уравнением теплопроводности 
для твердой фазы [7, 8]. 


Для приведения системы уравнений к безразмер- 
ному виду использовались следующие соотношения: 
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где х,у - координаты; и,ѵ - скорости по осям х,у\ 
Х,У и II, V - безразмерные координаты и скорости, 
соответствующие х,у и и,ѵ; I - длина полости по 
оси х\ - масштаб времени; г - безразмерное время; 
Ѵ іп — масштаб скорости (скорость входного пото- 
ка); Ѳ - безразмерная температура; Г 0 - температу- 
ра жидкости и твердого тела в начальный момент 
времени; Т іп — температура входного потока; у/ - 
функция тока и щ - ее масштаб; - безразмер- 
ный аналог у/; со - вихрь скорости и щ — его мас- 
штаб; О - безразмерный аналог со. 

Соответственно безразмерные уравнения 
Навье-Стокса в приближении Буссинеска в пере- 


115 




Известия Томского политехнического университета. 2006. Т. 309. № 5 


менных «вихрь скорости - функция тока - темпе 
ратура» для жидкой фазы (режим смешанной кон 
векции) и уравнение теплопроводности для твер 
дой фазы будут иметь вид: 
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здесь Сг = , — - число Ірасгофа; р - темпера- 

ѵ 

турный коэффициент объемного расширения; § - 
ускорение, создаваемое массовыми силами; 
2 уь 

Ке = - число Рейнольдса; Рг = — - число 

ѵ а 

Прандтля; Ро = ^- число Фурье; а - коэффици- 
ент температуропроводности. 

Начальные условия: 

У(Х,У, 0) = 0, 

П(Х,У,0) = 0, 

Ѳ(Х,У, 0) = 0. 

Граничные условия: 

- на внешнем контуре рассматриваемой области 
(кроме свободной поверхности жидкости) в од- 
ном варианте задаются граничные условия пер- 
вого рода Ѳ=0, а в другом - второго рода 

- на свободной поверхности жидкости: 
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на внутренних границах раздела твердой и жид- 
кой фаз, параллельных осям 07(0Я): 
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- на участке ввода жидкости: 
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- на участке оттока жидкости: 


дѲ(Х,Г) 


дУ 

й'Р 

07 


= 0, 


= 1 ; 


здесь Кі = число Кирпичева; Ві = — 

К(7 ~Т п ) Р Я 

- число Био; п - координата по нормали к рассма- 
триваемой поверхности; а - коэффициент тепло- 
обмена между внешней средой и рассматриваемой 
областью решения по границе у= 0; Т е - температу- 
ра окружающей среды; А„, Я/ - коэффициенты те- 
плопроводности твердой и жидкой фаз; Я ж/ - отно- 
сительный коэффициент теплопроводности. 

Задача (1-4) с соответствующими граничными 
и начальными условиями решена методом конеч- 
ных разностей [7]. 

Уравнения решались последовательно, каждый 
временной шаг начинался с вычисления поля тем- 
пературы в газовой и в твердой фазах, затем реша- 
лись уравнения Пуассона (3) для векторного по- 
тенциала. Далее определялись граничные условия 
для компонент вектора завихренности и решались 
уравнения движения. 

При решении системы уравнений (1-3) в каче- 
стве начального приближения использовались ста- 
ционарные решения, полученные в задачах при ма- 
лых числах Рейнольдса. 


4. Анализ результатов численного моделирования 

Численные исследования были проведены при 
следующих значениях: 

- безразмерных величин: Рг=7,1, 100<Ке<1000, Сг<10 5 ; 

- температур Т 0 = 293 К, 300<7)„<343 К, 273<Г<293 К. 

Диапазоны изменения безразмерных величин и 
размерных характеристик выбирались исходя из 
следующих условий. Первое - соответствие диапа- 
зонам изменения параметров в реальных вариантах 
реализации таких систем, как «жидкость - твердое 
тело - газообразная внешняя среда». Второе - ре- 
альные масштабы существования как вынужден- 
ной, так и естественной конвекции без вырожде- 
ния какого-либо из механизмов теплопереноса. 
Кроме того, контролировались условия геометри- 
ческого подобия реальных систем и используемой 
при математическом описании процесса модели. 
Также рассматривались не самые благоприятные 
для реализации «сопряженного» теплопереноса [8] 
диапазоны изменения как размерных величин, так 
и безразмерных комплексов. Так, например, увели- 
чение температуры вводимой в рассматриваемую 
полость жидкости, например, воды на 20.. .30° мо- 
гло существенно изменить (при прочих равных 
условиях) соотношение между характерными вре- 
менами реализации процесса вынужденной и есте- 
ственной конвекции, но экстремальные варианты 
не рассматривались. 
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